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ABSTRACT 

Context. The mass-loss process in Mira stars probably occurs in an asymmetric way where dust can form in inhomogeneous circum- 
stellar molecular clumps. Following asymmetries along the pulsation cycle can give us clues about these mass-loss processes. 

Aims. We imaged the Mira star X Hya and its environnement at different epochs to follow the evolution of the morphology in the 
continuum and in the molecular hands. 

Methods. We observed X Hya with AMBER in J-H-K at low resolution at two epochs. We modelled squared visibilities with geo¬ 
metrical and physical models. We also present imaging reconstruction results obtained with MiRA and based on the physical a priori 
images. 

Results. We report on the angular scale change of X Hya between the two epochs. ID CODEX profiles allowed us to understand 
and model the spectral variation of squared visibilities and constrain the stellar parameters. Reconstructed model-dependent images 
enabled us to reproduce closure phase signals and the azimuthal dependence of squared visibilities. They show evidence for material 
inhomogeneities located in the immediate environment of the star. 

Key words, stars: individual: HD 83048 - stars: imaging - stars: AGB and post-AGB - circumstellar matter - stars: mass-loss - 
techniques: high angular resolution 


1. Introduction 


When reaching the tip of the asymptotic giant branch (AGB), 
low-to-intermediate-mass stars begin a Mira phase before evolv¬ 
ing to the planetary nebulae and finally white dwarf stages. 
During this phase, their pulsation period of a few hundred days 
is associated wi th a several-magn itude photometric variability 
in the infrared dLe^ertrd, [19 93|) an d a mass-loss rate up to 
10 '^MQ/year dRamstedt et al.L 120081) . While this phenomenon 
is well explained by the A--mechanism, the way the radial oscil¬ 
lations interact with the stellar outer atmosphere and the way 
mass is dispersed into the interstellar medium (ISM) is not 
well understood. A molecular layer or MOLsphere extending to 
about 1 stellar radiu s from the photosphere has been systemati¬ 
cally observed (e. g. Mennesson et aTl 120021: iPerrin et '^l2004t 


iLe Bouauin et al.L 2009HWittkowski et al.Ll2068r 201 ll). The hy¬ 

pothe sis derived from spectroscopic H-band data (iHinkle et al.L 
1 19841) assumes free-falling material disassociated by the pulsa- 


* Based on observations obtained with the ESO VLTI/ATs telescopes 
under the program ID 084.D-0326. 


tion shock-front and recombining in the post-shock zone. For 
oxygen-rich stars like the one we study in this paper, one sce¬ 
nario predicts the formation of Fe-free silicate grains that allow 
a strong enough radiative pressure to drive a wind provided the 
grain size is about 1 micron , whic h is a typical size of grains 
observed in the ISM dHofneri [20081) . 


Another p iece of the M ira 
Ragland et al.l 


asymmetries. _ __ 

IWittkowski et al.l '( 2011 ) and ICruzalebes et al.l d2015] 


sta r 
(12005). 


puzzle are th e brigh t 
jLacouret^ d20()^ . 

I 


showed 

that asymmetries are common in Mira atmospheres and can be 
modelled (though not uniquely) as bright features containing 
up to several percent of the total flux. A number of phenomena 
have been proposed to explain these observations such as 
variable extinction resulting from condensation of dust grains 
in asymmetric regions, shock features or inhomogeneities in 
the molecular layers. Moreover, since protoplanetary nebulae 
(PPN) are a lso k nown to present asymmetric morphologies 
dSahai et al.L l2()()^ . we hope to understand the connection be¬ 
tween the AGB and PPN structures by studying the morphology 
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of the Mira circumstellar envelopes which are likely expelled 
with some degree of asymmetry. 

The goal of our study is therefore to characterize these asym¬ 
metries and establish a link with the mass-loss process. Near-IR 
spectro-interferometry is ideally suited for this study, since it de¬ 
livers spectral and spatial information on both the stellar pho¬ 
tosphere and the atmosphere or envelope. We here investigate 
the variation in the morphology of the Mira star X Hya based 
on two-epoch AMBER measurements. Previous interferometric 
obs ervations of X Hya w ere made with the IOTA interferome¬ 
ter (iRagland et al.L l2006l) in H band. They found that the object 
showed no sign of asymmetry up to a maximum baseline of 24.5 
meters for a visual pulsation phase of 0.5. Medium-resolution 
K-b and AMBER data were also obtained on X Hya at phase 
0.7 (IWittkowski et al.L 1201 Ih . Asymmetries were detected, more 
prominently in H 2 O bands, and were interpreted as clumps or 
inhomogeneities in this molecular layer. 

We here aim at imaging these asymmetries and observe their 
evolution along the pulsation cycle. In Sect. |2] we present the 
observations. In Sect. [3] we show the results of the parametric 
model htting analysis as well as radiative transfer modelling. 
In Sect. 0] we show the results of the image reconstruction ap¬ 
proach. Before concluding, we discuss these results with similar 
observations reported in the literature in Sect. |5] 

2. Observation and data reduction 

We observed the Mira star X Hya with AMBERh-EINITO 
(iPetrov et al.L l2007h at two epochs. In Pig. [1] these observation 
epochs are marked in the brightness variation cycle taken from 
the AAVSOQ database. They correspond to visual phases of 0.0 
and 0.2. X Hya is a M7e Mira star with a pulsat ion period of 301 
days, it is located at a distance of ~ 440 pc (IWhitelock et al.L 
I 2 OO 8 I 1 . A multi-wavelength observation is required to distinguish 
the molecular layer from the stellar photosphere. We thus per¬ 
formed a multi-wavelength J-H-K observation (covering wave¬ 
lengths from 1.1 yum to 2.4 yum) at low spectral resolution with 
AMBER. A log of the observations is presented in Table [T] The 
U-V plane we covered with our measurements was very similar 
at the two epochs as shown in Pig. |2] 

Reduced squared visibilitie s and closure phases were ob- 
tained with the amdlib package (iTatulli et al.L 1200'^ IChelli et al.L 
l2009l) . To optimise the stability of the transfer function, we ad¬ 
justed the frame selection with a threshold of 15 microns for the 
piston value. The percentage of the frames left after a selection 
based on the signal-to-noise ratio (S/N), and consequently the 
total integration time per calibrated point, varied for each run 
because of different atmospheric conditions. These information 
are provided in Table [T] 

The spectral calibration was based on a ht of telluric features 
and on the positions of flux minima between absorption bands 
(see Pig. Ell. We compared these positions between the observed 
spectra of the calibrators and a synthetic spectrum. This led us to 
multiply the wavelength table of our data by a factor of 1.04 to 
match the positions of the observed spectra with those of the syn¬ 
thetic spectrum. Interferometric calibrators were observed and 
their measurements bracketed those of X Hya. 

3. Model fitting of the squared visibilities 

The datasets we obtained at each epoch are rich and complex. 
Even at short baselines, squared visibilities (V^) show a signih- 

* http://www.aavso.org/ 


cant departure from a simple geometry such as an uniform disk 
model (Pig. |3ll. To understand the spectral and morphological 
signatures on the individually, we hrst divided the data into 
six spectral bands, as shown in Table|2] This division is based on 
an analysis of the visibility variation with the wavelength (see 
Pig. El which is coherent with a previous spectro-interferometric 
study of X Hya (IWittkowski et al.L 1201 ill and also based on the 
molec ular absorption bands shown in Pig. 4 of iLancon & Woo3 
(l2000f) . Although we labelled these bands ’’Continuum”, ”CO”, 
and ”7/2(9” based on the dominating ’’species”, they may also 
be contaminated by other molecular and atomic features that we 
cannot distinguish with a low spectral resolution. In this context, 
the label ’’Continuum” does not exclude some molecular con¬ 
tamination that may be present at these wavelengths, in particu¬ 
lar in the case of the Continuum (1) band, which overlaps with 
a part of an H 2 O absorption band. Our choice represents a com¬ 
promise between the location of molecular absorption bands and 
a relevant subdivision of the data to show their spectral variation. 

In the following, we present a progressive modelling ap¬ 
proach of the squared visibilities in the six spectral bands, which 
are representative of the centro-symmetric intensity distribution 
of the object. Before using a complex dynamical model, we per¬ 
formed a first standard modelling with the uniform disk function 
which delivers a basic reference size estimate. 


Spectral band 
number 

Wavelength range 

Contribution 

1 

1.1 - 1.4pm 

Continuum (1) 

2 

1.4- 1.5 pm 

HzOd) 

3 

1.5 - 1.7 pm 

CO(l) 

4 

1.7 - 2.1 pm 

H 2 O (2) 

5 

2.1 - 2.25 pm 

Continuum (2) 

6 

2.25 -2.4 pm 

CO (2) 


Table 2. Division of the covered spectral range in six spectral 
bands. 


3.1. Uniform disk model 

We hrst modelled a uniform disk (UD) to obtain a hrst estimate 
of the angular size of the object in the different bands and at the 
different position angles (PA) we observed. We only considered 
squared visibility data whose spatial frequencies belong to the 
hrst lobe so as not to bias this estimate with any contamination 
from limb darkening or structures smaller than the stellar diam¬ 
eter. An average of the measurements over the PAs is shown in 
Table[3] 

This hrst basic modelling description contains rich informa¬ 
tion. The UD diameter systematically increases between the two 
epochs in every spectral band. This agrees with the visual pulsa¬ 
tion cycle seen in Pig. 0] as the envelope expands, the tempera¬ 
ture decrea ses and corresponds to a dimmer visual magnitude. 

VINCI (iKervella et al.L l2()()(ih measurements of X Hya were 
reduced and analysed for comparison. They were taken in K 
band on 2002-01-27, corresponding to a phase of 0.25 which 
was close to the phase of the second epoch of our AMBER mea¬ 
surements (0.2). The UD diameter we derived is 7.55 + 0.70 
mas which is compatible with our results for band 5 (2.1-2.25 
jum) and 6 (2.25 -2.4 fim) with the overlap of error bars. 

The UD diameters reported in IWittkowski et al.l (1201 ih that 
were obtained for phase <I) = 0.7 (see last column of Table [3]) 
are also compatible with our measurements for bands 4 and 5 
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AAVSO DATA FOR X HYA • WWW.AAVSO.ORG 



Fig-1 - AAVSO V-band light curve of X Hya. The period of its brightness variations has been estimated to 301 days (IWhitelock et all 
l2008h . The two observation epochs are marked by red vertical lines and correspond to visual phases 0.0 and 0.2. 




u/A, (mas”'') 


Fig. 2. Uv coverage of the observations of the two epochs; left epoch A, right epoch B. Colours indicate the six spectral bands 
explicated in Table|2](l: red, 2: orange, 3: light green, 4; dark green, 5: turquoise, 6: blue). 


{H 2 O (2) and continuum). This is not the case for band 6 (CO 
(2)), however, for which Wittkowski et al. found a high value of 
9.0 + 0.1 mas. 


Linear limb-darkened disk diameters were also computed to 
m odel the data up to the second lobe (see Eq. 5 with B - 0 
in iHaubois et al.L l2009h . Results are presented in Table H) As 
judged by the high reduced values (whose maximum was 
100) and the extreme values found for the linear coefficient 
(forced to be between -1 and 1), the addition of limb darken¬ 
ing did not improve the modelling compared to the UD model. 
This suggests that data can only be reproduced by a more 
complex intensity distribution. 


3.2. Radiative transfer modeiiing 

To better understand the spectral dependence of this complex 
dataset, we used a grid of dynami c model atmosphere serie s 
computed with the CODEX code (llreland et al.L l2008l 1201 ih . 
These models develop a self-excited pulsation mechanism and 
sample opacities over 4300 wavelengths running from 200 nm 
to 50 fj.m, with 1 nm sampling between 200 nm and 3 /rm. Eor 
more details about CODEX and the grid of stellar p arameters for 
which the model series were computed, we refer to llreland et al.l 
(1201 ih . Eor each set of stellar parameters, these models pro¬ 
duce radial profiles of the intensity across the chosen spectral 
band that we turned into by a Eourier transform. These latter 
were finally fitted to our data by adjusting the Rosseland diame¬ 
ter and the pulsation phase. Eigure|4] shows V^-fitted curves for 
two triplets of baselines at two epochs. At the difference of the 
simple geometric models, these physical models that take into 
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Spectral 

band 

number 

(f>uD (mas) 

0 = 0.0 
(epoch A) 

(f>VD (mas) 

O = 0.2 
(epoch B) 

M<Pud)I 4>ud 

(puD (mas)' 

C) = 0.7 

1 

5.94 ±0.15 

7.41 ± 0.40 

0.25 ± 0.05 

- 

2 

6.90 ±0.13 

7.85 ± 0.49 

0.14 ±0.05 

- 

3 

6.03 ± 0.59 

7.19 ±0.85 

0.20 ±0.13 

- 

4 

6.52 ± 0.23 

8.24 ± 0.86 

0.27 ± 0.09 

8.3 ±0.1 

5 

6.20 ± 0.35 

7.02 ± 0.86 

0.14 ±0.10 

6.2 ± 0.1 

6 

6.80 ± 0.42 

7.54 ± 0.79 

0.11 ±0.09 

9.0 ±0.1 


Table 3. Fitting results of uniform disk models performed in 
the first lobe of the squared visibilities. Diameter values were 
averaged over the PAs. Results are expressed as yu + Icr, where 
yu is the average and cr is the standard deviation computed over 
all the PAs. The fourth column represents the relative variation 
of UD diameter between epochs A and B, rela tively to epoch A. 
'The l ast column reports measurements from IWittkowski et al.l 

dMIIl). 


account H 2 O and CO opacities do reproduce the variation of 
with wavelength fairly well. The results of the model parameter 
fitting are summarised in Table |5] 


Models 

Epoch A, O = 0.0 

Epoch B, O = 0.2 

Best R52 model 

388040 

388040 

dposs (mas) 

4.02 

4.93 


8.1 

21.9 

^CODEX 

0.7 

0.7 

Best o54 model 

262160 

286060 

SRoss(mas) 

4.76 

6.03 

2 

red 

lA 

19. 

^CODEX 

0.7 

0.4 


Table 5. Results of the fitting of the CODEX radiative transfer 
models. 


As for the previous parametric modelling, the Rosseland di¬ 
ameters found with the model fitting all increase from epoch A 
to epoch B. Relatively to epoch A, the increase in diameter from 
epoch A to epoch B was 23 % and 27 % for the R52 and o54 
CODEX series, respectively. Com pared to the same mode lling 
for O = 0.7 reported in Table 3 of IWittkowski et all (1201 ih . we 
globally find a poorer fit by the attempted models to the data as 
can be seen from the values of the reduced The R52 model 
fits the data almost equally well as the o54 series. These two 
models differ from each other in the parent star parameters (o54: 
M = 1.1 Mo, L = 5400 Lo, R = 216 Rq, P = 330 days; R52: M 
= 1.1 Mo, L=5200 Lo, R= 209 R©, P= 307 days). The fact that 
the CODEX series were computed for a limited parameter space 
could also explain the lack of a better fit to our data. Moreover, 
the phases fitted by the CODEX models (0.7 and 0.4) do not 
match the visual phases (0.0 and 0.2, respectively). This discrep¬ 
ancy in the phases might be explained by the degeneracy in the 
CODEX models with the parameters and the pulsation phases; 
the same atmosphere might be equally modeled by two differ¬ 
ent sets of parameters and phases. On the other hand, CODEX 
models are ID and do not take into account non-radial features. 
The presence of such non-radial features in some azimuth ranges 
causes the visibility curve to depart from a pure radial profile 
and would explain why the phase determination is not accurate. 
Among the possible origins for non-radial intensity distributions 
we can cite : surface spots, non-radial pulsation, an asymmet¬ 
ric shock front, an inhomogeneous molecular layer and dust en¬ 


vironment, presence of a companion, etc. In the following, we 
present an imaging work that aims at unveiling the non-radial 
and asymmetric properties in the morphology of the object. 

4. Imaging 

In the previous sections, we only took into account the data, 
which means that the object was thought to be centro-symmetric. 
However, since AMBER recombines the beams from three tele¬ 
scopes, we do have a phase information called closure phase that 
is a measurement of the asymmetry in the brightness distribution 
of the object. Eigure |5] shows that closure phase departs signifi¬ 
cantly from the 0 modulo 180° values which means that X Hya 
was not centro-symmetric in most of the spectral bands. 

The closure phase signal of epochs A and B exhibits a noisy 
pattern at the spectral band 1, which is located around 100m of 
maximum projected baseline (second visibility lobe). This might 
be due to a combination of poor atmospheric conditions and low 
object visibilities at these spatial frequencies. The object can 
be considered centro-symmetric below ~80 meters which cor¬ 
responds to ~5.5 mas at 2.15 fim. Eor epoch B, closure phases 
show a much higher dispersion and higher uncertainties. The 
maximum baseline where the closure phases depart from the 
null value lies between 55 and 60 meters, which indicates a 
higher degree of asymmetry than for epoch A. These results 
agree w ith the //-band measurements presented in iRagland et al.l 
(l2006h that reported that the closure phase signal of X Hya was 
consistent with the null value up to a maximum baseline of 24.5 
meters (observations at d) = 0.5). 

To reconstruct the spatial intensity distribution out of such 
a complex dataset of and closure phases, it is ma ndatory to 
use im age reconstruction algorithms such as MiRA dThiebaud 
l2008h . However, given our limited uv coverage in individual 
bands, we cannot aim at reconstructing a model-independent im¬ 
age. Nevertheless, we carried out an image reconstruction ap¬ 
proach based on our parametric modelling results. We therefore 
used ID intensity profiles from the CODEX models that we 
turned into initial guess and a priori image for the reconstruc¬ 
tion. Several regularisation techniques were used including the 
total variation (TV) regularisation recommended i n benchmark¬ 
ing st udies for image reconstruction with MiRA (iRenard et all 
I 2 OI Ih . The best results were obtained with the TV and quadratic 
regularisations towards an a priori object image derived from 
the R52/388040 model intensity profile (Eig.|6]l. 

Eor the quadratic regularisation, the reconstructed images 
were obtained by adjusting the data and by regularising towards 
this prior image. This was done by minimising the following cri¬ 
terion: 

^ l^fpriori^^y)] ( 1 ) 

-c.y 

with 

/prior = ^ [I{x, y) - P{x, y)f (2) 

x,y 

where I(x,y) is the reconstructed image and P(x,y) is the prior 
image, x and y are the coordinates, and p, also called hyperpa¬ 
rameter, sets the weight between data (V^ and closure phase) 
fitting, represented by the criterion Jdata, and regularisation to¬ 
wards the prior image, the criterion called /prior- 

Eor the TV regularisation, the a priori object image was 
used as a guess image to guide the data fitting. As shown in 
Eigs.|7l[8]|9][l0l[II] and[T2l data are convincingly fitted after the 
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image reconstruction process. Table |6] summarises the reduced 
values as well as hyperparameter values and functions used 
in MiRA to reconstruct images. 


Spectral band 
number 

Regularisation 

C) = 0.0 
(epoch A) 

Regularisation 

0 = 0.2 

(epoch B) 

1 

Quadratic / le5 / 3.6e-2 

Quadratic / le6 / 0.3 

2 

Quadratic / le3 / 1.2e-3 

Quadratic / le5 / 4.2e-2 

3 

TV / le2 / 0.25 

Quadratic / le5 / 1.8e-2 

4 

Quadratic / le4 / 1.2e-2 

TV / le3/ 1.4 

5 

Quadratic / 5e3 / le-3 

Quadratic / le3 / 1.1 

6 

TV / le2 / 0.4 

TV / le2/0.5 


Table 6. Regularisation functions and hyperparameters used in 
MiRA to produce reconstructed images presented in Figs. [T3lfT4l 
and [15] Reduced values are specified and represent the value 
of the minimised criterion per datum at the end of the recon¬ 
struction process. 


Reconstructed images for the six spectral bands at the two 
epochs are presented in Figs. [T3lfT4l and[T5] 

The fact that the uv coverages at the two observing periods 
are almost identical suggests that the morphology of X Hya re¬ 
ally changed from epoch A to epoch B. At a 25%-intensity level 
of the reconstructed images shown in Fig. [13] to [15] the stellar 
morphology is disk-like in all spectral bands. If we use this in¬ 
tensity threshold to estimate the stellar diameter, the relative in¬ 
crease of this value between epochs A and B ranges from 13% to 
18% depending on the spectral band (see Table [Tj. Considering 
the overlap of 1 cr-error bars, these results agree with the re¬ 
sults of the UD modelling shown in Table [3] except for band 1 
for which the 25%-intensity diameter variation is found to be 
weaker than the UD diameter variation. A difference between 
the two diameter estimator ( 025 % and ipun) values is expected 
because they take into account different intensity levels for the 
same given image. However, the relative variation of the two es¬ 
timators should give close values if most of the spatial intensity 
distribution retains the same profile between the two epochs. For 
band 1, it seems to indicate that this is not the case. This could 
either be a real astrophysical effect or a consequence of noisy 
data (particularly for epoch B, right part of Fig. [7]i in the recon¬ 
structed images (Fig.HjI). The latter case could be explained by a 
combination of low visibilities and the effect of atmospheric tur¬ 
bulence, stronger at these wavelengths (7 band), on the accuracy 
of both and closure phase data. 


Spectral 

band 

number 

025 % (mas) 

0 = 0.0 

(epoch A) 

025 % (mas) 

0 = 0.2 

(epoch B) 

A(025%)/025% 

1 

5.01 ±0.15 

5.67 ±0.19 

0.13 ± 0.04 

2 

4.95 ±0.16 

5.64 ±0.17 

0.14 ± 0.03 

3 

4.76 ± 0.21 

5.57 ± 0.23 

0.17 ± 0.05 

4 

4.73 ± 0.20 

5.57 ± 0.23 

0.18 ±0.05 

5 

4.84 ±0.21 

5.57 ± 0.23 

0.15 ±0.05 

6 

4.90 ±0.18 

5.57 ± 0.23 

0.14 ± 0.04 


Table 7. Stellar diameters estimated from the 25% iso-intensity 
levels of Fig. [T3lfT4l and[T5] Diameter values were averaged over 
the PAs. Results are expressed as yu + Icr, where is the average 
and cr is the standard deviation computed over all the PAs. The 
fourth column represents the relative variation of stellar diameter 
between epochs A and B, relatively to epoch A. 


For the two epochs, the shape of the star varies quite sig¬ 
nificantly between pairs of conjugated spectral bands (1-5, 2- 
4, and 3-6). The first reason is the fact that each spectral band 
corresponds to a different number of spectral channels and 
consequently to a different number of and closure phase 
points (asymmetry information) for the image reconstruction. 
Secondly, the size of the interferometric beam, which is larger 
at longer wavelengths, makes it more difficult to resolve small 
features (including asymmetries) at longer wavelengths. Finally, 
as can already be seen from Fig. [4] the size of the object as pre¬ 
dicted by the radiative transfer calculations, does vary for all the 
bands. This explains why we can see different morphologies be¬ 
tween bands 1 and 5 (continuum bands), 2 and 4 (H 2 O) and 3 
and 6 (CO) for each epoch. 

We can distinguish two main components in the recon¬ 
structed images (Figs.[T3]to[T5]i: 

- the stellar disk, whose diameter variation is in agreement 
with previous estimations of geometrical models, and whose 
intensity is above 25% of the maximum intensity level. 
Above this level, all images are centro-symmetric apart from 
the Continuum (1) at Epoch B see the discussion above. 

- the close environment or atmosphere, whose intensity is be¬ 
low 25% of the maximum intensity, that shows heteroge¬ 
neous morphologies and whose characteristic angular size 
can reach up to 10 mas. 

We add a word of caution at t his point. The 201 4 
Interferometric Imaging Beauty Contest (iMonnier et al.L l2014l) . 
which aimed at comparing image reconstruction methods, 
showed that reconstructed images based on the same dataset di¬ 
verge for intensities lower than 25% of the maximum relative 
intensity (see Fig. 5 of the quoted paper). Below this intensity 
level, it is therefore difficult to distinguish intrinsic astrophysi¬ 
cal features from image reconstruction artefacts when one uses 
a single algorithm. Nevertheless in the following, we speculate 
on the possible nature of features and asymmetries if they prove 
to be real astrophysical properties of the star’s morphology. 

This hypothesis is supported by the fact that for most of 
the reconstructed images presented here, the two regularisation 
types and various hyperparameter values revealed the same fea¬ 
tures for intensities above 10% of the maximum intensity (// val¬ 
ues typically ranged from le3 to le6 for the quadratic regularisa¬ 
tion and from lei to le3 for the TV regularisation). Exceptions 
are the CO(2) images at epoch A and the Cont(l) and H20(2) 
images at epoch B, which exhibit different environment shapes 
for the two regularisation methods even though they all pro¬ 
duce good data fitting, such as the one presented in Table [6] 
Nevertheless, for all spectral bands, the maximum size of the 
environment and the shape of the stellar disk remain stable with 
different regularisation methods. 

Eor the two epochs, Eig.[T6]shows intensity contour lines of 
the CO and H 2 O environments overlaid on the continuum im¬ 
ages. Upper panels correspond to a smaller dirty beam and are 
better indicated to study the extension and the shape of the envi¬ 
ronment. Eor the two epochs, the azimuth-averaged extension of 
the environment is smaller in the continuum than in the H 20 (l) 
band, the environment seen at CO(l) being included between 
the two. However, due to the irregular shape of these structures, 
these relative extensions vary for given azimuths. 

The impact of the uv coverage (represented as the dirty 
beam in the lower right corner of each image) on the environ¬ 
ment’s shape is well visible on some images. Particularly for 
Continuum(2) (epoch B) and CO(2) (epochs A and B) and for 
H20(2) at epoch B, we can see a global SE-NW elongation that 
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is almost aligned with the dirty beam (which is larger at these 
wavelengths). If this elongation has no any obvious physical ori¬ 
gin for a pulsating star like X Hya, it remains uncertain, however, 
whether it is entirely due to an artefact of the reconstruction. 

This is not obvious for H20(2) at epoch A where the 
curve (Fig. [Toll shows a clear departure from a simple radial sym¬ 
metry. Moreover, the elongation is not exactly aligned with the 
dirty beam which potentially means that it is not an artefact. 

We now focus on the CO(l) and H20(l) bands that show the 
most significant closure phase and irregular Y^ data: 

- H 20 ( 1 ): although the dataset contains more points and has 
larger error bars for Epoch B, it looks more complex than 
at Epoch A as can be seen from the non-regular curve 
and the scattering in the closure phase signal. Comparing the 
images at the two epochs, one can see the general morphol¬ 
ogy remained similar. It seems that the stellar surface and 
the close environment of the star expanded from Epoch A to 
Epoch B, and that some asymmetries appeared in the close 
environment at Epoch B (Eig.[T4li. 

- CO(l): on the two images, the close environment doesn’t 
present the same structure. If the elongation of Epoch A is 
probably an effect of the uv coverage, it doesn’t appear at 
Epoch B. Instead, two arms reproduce correctly the two parts 
of the second lobe. 


5. Discussion 

The main goal of this observing campaign is to follow the evolu¬ 
tion of asymmetries during the pulsation cycle of a Mira star. The 
degree of asymmetry as seen by the values of the closure phase 
signal increased from epoch A (maximum of the visual pulsa¬ 
tion cycle) to epoch B (see Eig. |5]l, where the system star-H envi¬ 
ronment is also generally more extended angularly. We mention 
here that instead of the closure phase, an alt ernative estimator 
called the centrosymmetry parameter (CSP, ICruzalebes et al.L 
I 2 OI 4 I) would be also interesting to use since it is designed to be 
more sensitive to asymmetries than the closure phase. 

Compared to a sim ilar imaging campaign of T Lep 
(iLe Bouauin et all l2009l) . it is noteworthy that we did not find 
any discontinuity in the intensity distribution of X Hya between 
the stellar surface to the close molecular environment. Moreover, 
CODEX models do not predict such a discontinuity in general. 

The asymmetries in the close environment (typically 1.5 stel¬ 
lar radius) in H 2 O and potentially CO molecular bands can be 
explained by an inhomogeneous distribution of circumstellar 
matter. We can therefore support the interpretation of previous 
medium-resolution AM BER observations of X Hya reported in 
IWittkowski et al.l (1201 ill , who suggested the presence of clumps 
in the water vapor band. Extrapolating from the H20( 1) image at 
epoch B (Eig.[T4l). clumps would have a typical angular scale of 
~ 2 mas, representing a few percent of the total flux each. This 
conclusion agrees with the analysis of the c losure phase signal of 
R Cnc reported in IWittkowski et akl (1201 ll) . which also suggests 
the presence of unresolved spot (up to ~ 3 mas) contributing to 
up to ~ 3% of the total flux in the H20(2) band. It is noteworthy 
that this ~ 2-3 mas clump size also represents the typical thick¬ 
ness of the atmosphere as modelled by CODEX series where the 
typical features of the stellar spectrum are formed. 

Clumps in the atmosphere are indeed more likely to be 
present after the visual maximum phase, that is after the shock 
wave had time to penetrate the most external layers and pos¬ 
sibly enhanced locally inhomogeneous distributions of matter. 


On the other hand, recent progress in 3D radiation hydrodynam¬ 
ics simulations revealed that convection alone can account for 
asymmetrically extended features (synthetic images in J , H an d 
K bands are presented in Eig. 5 of IChiavassa & Erevtael 120141) . 
It therefore remains possible that asymmetries also have a con¬ 
vective origin, whether they are formed directly by convective 
motions or whether convection set up density inhomogeneities 
and velocity fields on the scale of the photospheric radius, and 
these could be propagated into the environment by pulsation. A 
detailed comparison of intensity levels of the asymmetric fea¬ 
tures between these simulations and our multi-spectral observa¬ 
tions has to be carried out in a future work. 

6. Conclusions 

We reported on AMBER low-resolution observations of the 
Mira star X Hya at two epochs. Using advanced radiative trans¬ 
fer modelling tailored to studying dynamics in Mira star at¬ 
mospheres, we successfully modelled the spectral variation of 
squared visibilities. Model-dependent images show asymmetries 
that are located in the immediate environment (at ~ 1.5 stellar ra¬ 
dius). As predicted by the CODEX models, no discontinuity is 
observed between the surface and the immediate environment 
of the star. The two best-fitting models provided an estimate of 
the fundamental parameters of X Hya: M = 1.1 M© (the mass 
is the same for the two models), L = 5300 + lOOL©, R = 212 + 
3 Rg and P = 319+12 days. More epochs are needed to prop¬ 
erly model the atmospheric dynamics and test further radiative 
transfer codes. More parameter spa ce should be explore d for the 
CODEX series as also concluded in iHillen et al.l (l2012h . 

We resolved asymmetric features on the Mira star X Hya: 

1. As shown by the evolution of the closure phase signal (see 
Sect. 4) the object intensity distribution is more asymmetric 
at epoch B (O = 0.2) than at epoch A (O = 0) 

2. Asymmetries may explain the discrepancy between the vi¬ 
sual phase and the phase derived from the (ID) centro- 
symmetric dynamical models of pulsating atmospheres we 
used. 

3. The asymmetries located in the immediate environment 
could be inhomogeneous material enhanced by the shock 
wave passage or a feature of a convective pattern. 

Eurther observing campaigns are necessary to confirm these 
conclusions on X Hya and on other Mira stars. A better uv 
coverage and more accurate closure phase measurements are 
needed to reconstruct model-independent images and unam¬ 
biguously estimate the characteristics of asymmetries. Ideally, 
these new datasets should also make use of various imaging 
reconstruction algorithms to efficiently distinguish reconstruc¬ 
tion artefacts from astrophysical asymmetries in the close envi¬ 
ronment. The hypothesis that asymmetries are due to a convec¬ 
tive pattern could then be efficiently tested by comparing multi- 
band images with 3D radiat ion hydrodynamics simulations (e.g. 
IChiavassa & Erevt^l2014li . 
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Date/MJD 

Configuration 

Calibrator 

Coherence time’ (ms) 

Fried parameter’ (cm) 

Erame selection 

Integration time (s) 
calibrated point 

Epoch A (Visual Phase: 0.0) 

09-01-2010/55205 

DO-HO-KO 

HR3802 

5.3 ± 1.2 

2.0 ±0.4 

20% 

48.5 

10-01-2010/55206 

DO-Gl-HO 

HR3802 

3.2 ± 0.8 

2.2 ± 0.5 

40% 

97.0 

Epoch B (Visual Phase: 0.2) 

16-03-2010/55271 

DO-HO-KO 

HR3802 

16.0 ± 6.5 

3.0 ± 1.2 

35% 

85.0 


DO-Gl-HO 

HR3749 

16.0 ± 6.5 

3.0 ± 1.2 

35% 

85.0 


Table 1. Log of the observations. ' As an indication on the atmospheric conditions, coherence times and Fried parameters were 
taken from the DIMM measurements in the visible and are expressed as /r + Icr, where /r is the average and cr is the standard 
deviation computed over each half-night of our observations. 



Projected baseline [m] Projected baseline [m] 

Fig. 3. Squared visibilities data at epoch A (left) and epoch B (right). Colours indicate the six different spectral bands as in Fig.|2 
The solid line represents a 6 mas (epoch A) and a 7 mas (epoch B) uniform disks at 1.75 micron (average wavelength of the covered 
spectral intervals). 


Spectral band 
number 

<I>LD 

0) = 0.0 
(Epoch A) 

Ald 

0 = 0.0 
(Epoch A) 

^LD 

O = 0.2 
(Epoch B) 

Ald 

C) = 0.2 
(Epoch B) 

1 

5.73 ± 0.04 mas 

1.00 ±0.11 

8.47 ± 0.08 mas 

1.00 ±0.12 

2 

8.56 ± 0.28 mas 

1.00 ±0.12 

6.88 ± 0.07 mas 

1.00 ± 0.10 

3 

5.67 ± 0.05 mas 

1.00 ±0.11 

6.96 ± 0.02 mas 

1.00 ±0.05 

4 

6.37 ± 0.34 mas 

1.00 ±0.11 

10.58 ± 0.10 mas 

0.58 ± 0.06 

5 

6.00 ± 0.37 mas 

1.00 ±0.10 

7.19 ± 0.06 mas 

-1.00 ±0.10 

6 

6.46 ± 0.36 mas 

-1.0 ± 1.16 

7.58 ± 0.05 mas 

1.00 ± 0.12 


Table 4. Fitting results of linear limb-darkened disk models. These measurements were made over all PAs. The LD coefficient was 
forced to be between -1 and 1. Reduced chi2 values are between 2 and 100. 
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Fig. 4. Modeling of the squared visibility data (in black) with R5 2 (red) and o54 (blue) C ODEX models. Our division in 6 spectral 
bands is indicated as well as the molecular absorption bands from iLancon & WoodI (l2000ll . Each panel shows three visibility spectra 
corresponding to the three baselines of one triplet: DO-HO, HO-KO, and DO-KO (upper panels) and DO-Gl, DO-HO, and Gl-HO 
(lower panels). The baseline ground lengths for DO-HO/HO-KO/DO-KO and DO-Gl/DO-HO/Gl-HO are respectively 64m/32m/96m 
and 72m/64m/72m. The left panels represent Epoch A whereas the right ones represent Epoch B. The position angle variation of 
the configuration between the two epochs is less than 1 degree for DO-HO-KO and less than 5 degrees for DO-Gl-HO. 




150 


Maximum projected baseline [m] 


Maximum projected baseline [m] 


Fig. 5. Closure phases in degrees at epochs A (left) and B (right). Colours indicate the six different spectral bands with the same 
coding as in Fig.|2] 
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Fig. 6. Image built from the CODEX intensity profile R52/388040 at 1.9 micron. This image was used as a prior and initial guess 
for the image reconstruction of band 4 datasets (H20(2)). 



Maximum projected baseline [m] 



Fig. 7. Squared visibilities (Upper panels) and closure phases (lower panels) data points and their fit by the image reconstruction 
process are plotted in red and green respectively, for Epoch A (left) and Epoch B (right). This dataset corresponds to band 1, i.e. 
Continuum (1). 
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Maximum projected baseline [m] 

Fig. 8 . Same as Fig. |7]for band 5, Continuum (2). 



Maximum projected baseline [m] 



Maximum projected baseline [m] 

Fig.9. Same as Fig.|7]for band 2, H 2 O (1). 
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Fig. 10. Same as Fig. [Tjfor band 4, H 2 O (2). 
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Fig. 11. Same as Fig. [Tjfor band 3, CO (1). 
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Fig. 12. Same as Fig. [Tjfor band 6, CO (2). 
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Fig. 13. Reconstructed images with MiRA for spectral bands 1 and 5 corresponding to the continuum bands at epochs A (left 
column) and B (right column). The interferometric beam size is displayed in the bottom right corner of each image. White contours 
represent iso-intensity levels of 25%, 75% and 95% of the maximum intensity of the image. 
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Fig. 14. Same as Fig.[T3for the H 2 O bands (bands 2 and 4). 
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Fig. 15. Same as Fig.[T3for the CO bands (bands 3 and 6). 
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Fig. 16. Continuum images overplotted (Fig.[T3ll with H 2 O (blue) and CO (white) intensity contours (from Figs.[T4landfT5ll. Intensity 
contour levels represent 2% and 5% of the maximum intensity of each image and allow comparing the extension and shape of the 
environment between spectral bands for epochs A (left panels) and B (right epanels). The upper panels represent the Cont (1), 
H20(1) and CO(l) bands, whereas the lower panels represent the Cont (2), H20(2) and CO(2) bands. 
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